New CLEO Results for \V cb \ and \V ub \ 

Roy A. Briere, representing the CLEO Collaboration 

Carnegie Mellon University; 5000 Forbes Ave.; Pittsburgh, PA 15213 

Abstract. We report recent measurements from CLEO of the first two moments of the photon 
energy spectrum for b — » sy decays and the hadronic recoil mass in B — > X c £v. These physical 
quantities allow one to fix non-perturbative parameters occurring in calculations based on HQET 
and QCD. Predictions for semileptonic decay rates within this same framework depend in addition 
on the CKM matrix elements V qq i governing quark mixing. We can thus extract \V c h\ from the 
inclusive semileptonic decay rate of B mesons, and \V„b\ from the lepton endpoint spectrum of 
B — > X u iv. Model dependence is reduced except for the assumption of quark-hadron duality. Finally, 
we update the classic measurement of \V c b\ from B — > D*£v at zero recoil. 

INTRODUCTION 

Most of the ad-hoc parameters of the Standard Model are contained in the flavor sector. 
Mysteries of mixing, mass generation, and CP violation all meet here. B physics of- 
fers the possibility of direct measurements of two of the Cabibbo-Kobayashi-Maskawa 
(CKM) mixing-matrix elements, \V c b\ and \V u b\, as well as access to two others, \V t d\ and 
\V ts \, via loop processes. We concentrate here on determinations of the former pair via 
semileptonic B decays. 

Experimentally, we measure the number of certain semileptonic decays; we can con- 
vert this to a branching ratio via knowledge of N B g (the number of B pairs present in our 
data sample) and finally to a partial width by using B lifetime measurements from else- 
where. Both these partial widths and certain other kinematic 'moments' may be calcu- 
lated in a systematic Heavy Quark Effective Theory (HQET) [[[], |2|] and QCD expansion. 
The expressions depend on some a priori unknown non-perturbative parameters, most 
notably A, A4 and X2. The kinematic moments will allow us to independently determine 
these parameters for use in a self-consistent way in other formulae. 

Expressions for semileptonic decay rates also depend on CKM matrix elements. 
One can thus extract \V c b\ from the inclusive semileptonic rate, or perform a more 
intricate extraction of \V u b\ from the rate at the lepton endpoint ofb-^u semileptonic 
decays. The partial width calculations rely on quark-hadron duality [|3p. That is, they 
assume that for sufficiently inclusive quantities, quark-gluon calculations can be used 
for observed hadronic processes. The chief issues are non-perturbative effects and the 
need to average over enough hadronic states, which is even more problematic for \V u b\ 
from the endpoint. It is therefore desirable to compare these results to those obtained 
from previous methods; one such result is also updated here. This is our determination 
of \V c b\ from B — > D*£v at zero-recoil as favored by more familiar HQET treatments. 

The CLEO II detector [|]] and the CLEO II. V upgrade with a silicon vertex detector 
[^]] and new drift chamber gas [^] are described elsewhere. Overall, in CLEO II and II. V 
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FIGURE 1. Left: Observed photon spectrum shown in a) with the scaled continuum background 
prediction from off-resonance data, and after continuum subtraction in b), where the BB background 
is now displayed. Right: Measured photon spectrum for b — » sy events. 



data, the T(45) to continuum luminosity ratio is about 2.1 and the effective BB cross 
section is 1.06 nb. The b — > sy and the B — > X u £v analyses both use the entire CLEO 
II and II.V datasets with a total luminosity of about 13.5 fb -1 and about 9.7 x 10 6 BB 
pairs, while the B — > X c £v and the B — > D*£\ analyses use only the CLEO II portion with 
about 3.3 x 10 6 BB pairs. 



THE b -> PHOTON SPECTRUM 

This final state is distinguished by a high-energy photon. Our new analysis uses a lower 
cut of 2.0 GeV on Ey, reducing model-dependence systematics. Continuum background 
dominates, but a factor of one hundred suppression is achieved by using event shape in- 
formation, kinematics of detected leptons, and a pseudo-reconstruction technique (seek- 
ing the best K(nn) y combination). A neural network combines information to give a sig- 
nal weight for each event. Both the raw photon spectrum indicating background sources 
and the final subtracted spectrum are displayed in Figure [l|. Though BB background 
increases at the lowest energies, we still see sensible behavior albeit with larger errors. 

Our focus here is extracting moments of the Ey spectrum; for other analysis details, see 
[0]. This spectrum, naively a sharp line, has a width determined largely by the b quark 
Fermi motion and somewhat by the varying recoil mass (i.e., QCD effects). Two smaller 
sources of width are the small known B boost (P ~ 0.06) and resolution smearing. 

We determine moments directly from the data, accounting for energy-dependent 
efficiency, resolution, and boost smearing. As a check, we also take the Ali-Greub 
[|]] or Kagan-Neubert models [Bp, propagated through our full GEANT-based detector 
simulation, and take moments of these models with parameters which best fit the data. 
All extractions are consistent, and variations are reflected in the systematic error. 



We find, for Ey > 2.0 GeV [|TOp 



<£ Y > = (2.346 ±0.032 ±0.011) GeV 
<£ 2 > - <E y > 2 = (0.0226 ±0.0066 ±0.0020) GeV 2 

where the brackets < ... > denote the average value. 

The first moment, again with Ey > 2.0 GeV, is calculated as [ [TT], |T2"| , |T3| ]: 
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+ 0(l/AfS)]. 

To lowest order, (Ey) = A [Mg — A] . The parameter A measures the energy of the 
light degrees of freedom: the 'brown muck' surrounding the heavy-quark. There are 
further parameters appearing in general at second order in 1 /Mg\ they are absent in this 
particular case, and are discussed later. Finally, the third-order parameters % and p ; are 
estimated as O (0.5 GeV 3 ) with variations included as systematics. 

From the first moment alone, we find that A = (0.35 ± 0.08 ± 0. 10) GeV [110]. Here, 



the entire experimental error is given first and the second error is due to the theoretical 
extraction. We avoid using second moments which give poorer determinations of the 
parameters and, in the case of hadronic moments below, have expressions which do not 
appear to converge as rapidly. Since teh value of A is not meaningful out of context, 
one must do all calculations consistently with respect to scheme and order. We choose 
MS, O ( 1 /Ml), O (Poa 2 ) everywhere. 



HADRONIC MOMENTS IN B -> X c £v 

This analysis uses both e and p with 1.5 < pi < 2.5 GeV/c and we reconstruct the 
neutrino properties via four-momentum balance with techniques developed at CLEO 



[|14J. Great care is taken to account exactly once for all observed particles; cuts on 
charge balance, a multiple-lepton veto, and {E^ iss — p^iss) help ensure that the missing 
four-momenta is due to a single missing neutrino. The resolution on missing momentum 
is o(pmiss) — HOMeV/c. After checking consistency with zero missing mass, we use 
E\ = \p m iss\ rather than E m i ss since it has better resolution. We use B decay kinematics 
to determine M| = A/| + M 2 ^ — 2EgE^ + 2pg ■ p^ without needing to observe the 
hadrons or group them into those from the B vs. the B decay. The small final dot 
product, which averages zero, is ignored. Thus, four-momentum balance of the entire BB 
event measures the neutrino properties, leaving the hadronic recoil system unobserved. 
A typical exclusive mode analysis would instead observe the hadron(s) and have an 
unobserved neutrino. We expect 95% b — ► civ after continuum subtraction; the rest 
(c — > slv secondaries, b — > uiv) is subtracted via Monte Carlo simulation. 



2000 



1850801-005 




1850701-004 



2 

M^(GeV 2 ) 



~i 1 r 

■Experimental 

□Total 




1.0 



FIGURE 2. Left: Observed recoil mass in B — » X c £v showing data as points along with a fit to D, 
D* and heavier charm meson contributions. Right: Constraints from measured b — > sy photon energy 
and B — » X c i\ recoil mass first moments in the A — A-i plane. The ellipse indicates A% 2 = 1, including 
systematic errors. 



Final results for the moments are calculated from the M\ distributions corresponding 
to the mixture of D£\, D*£v, and Xh£v spectra which best fit the data. We take moments 
of the generated M§ distributions while fitting the data to reconstructed quantities passed 
through the full physics and detector simulation and analysis, hence accounting for the 
B boost, resolution and efficiency. Heavy states Xh beyond the D and D* include D** 
states modeled with ISGW2 [ |T3| ] and non-resonant flWji treated with the Goity-Roberts 
prescription [ |16| ] . Their normalization is fixed by data. The fit distributions are shown in 
the left panel of Figure 0. 



We finally arrive at [|17Q: 

<M\-M 2 D > = (0.251 ± 0.023 ± 0.062) GeV 2 
<{M\-M 1 D ) 1 > = (0.639 ±0.056 ±0.178) GeV 4 
<(M|-<M|>) 2 > = (0.576 ±0.048 ±0.163) GeV 4 

where Md denotes the spin-averaged D,D* mass. The main systematics include the 
neutrino reconstruction efficiency and the models of theX# states. 



The theoretical expressions for the moments QTS] , [19], gOp use a consistent scheme and 
include (most of) the effects of the lepton energy cut. Unlike the mean photon energy 
discussed earlier, the second order HQET expansion parameters appear here. These are 
Xi, related to the Fermi motion energy of the b quark, and X2, measuring the QCD 
hyperfine splitting; the latter is fixed from mg* — mg as measured by others. The third- 
order terms are treated as before. 

Combining with the b — ► sy result, we find [[17]]: 



A = (0.35 ±0.07 ±0.10) GeV 

Xi = (-0.236±0.071 ±0.078) GeV 2 



The errors have the same meaning as for the b — > sy moments. The results are best 
viewed in the A — X\ plane; see the right panel of Figure 



THE INCLUSIVE SEMILEPTONIC RATE AND \V cb \ 

It is of course also possible to calculate the zeroth moment for semileptonic decays; this 
is simply T s i = T(b — > c£v) Pffi [2"3"| ]. The expression looks like a free-quark decay, 
akin to the classic muon decay rate, with a phase-space factor for finite m c , augmented 
by QCD corrections and the HQET expansion. The actual formula used, gleaned from 
calculations in © j24j, H, 0, 01, may be found in [ [T7| ] . 



Having consistently determined the lower-order HQET parameters with moments, we 
are poised to extract \V c b\. The inclusive semileptonic rate is taken from a venerable 



CLEO result [ J28| ] using the tagged di-lepton method []29|]. After subtracting 1% from the 
published value to correct for B -> X M £v, we have: S (fl — > X c iv) = ( 10.39 ± 0.46) %. 

We convert to r,/ using x s ± = (1.548 ±0.032) ps [[p, x B o = (1.653 ±0.028) ps JHJ], 
and /+-//00 = 1.04 ±0.08 g. Combining, we finally arrive at [0: 

IV^I = (4.04 ± 0.09 ± 0.05 ± 0.08) x 10" 2 

The largest errors are from (in order) the measurement of T s j, the HQET parameters 
A, A4, and the scale for Oy. This yields a precise (3.2%) determination, but given the 
global quark-hadron duality issues one should compare to results from other methods. 



EXTRACTING \V ub \ FROM THE B -> X M £v ENDPOINT 

One might ask, why not do another expansion in A, Xi , X2 for the B — > X M £v partial 
width? One can do this for the fully inclusive width; for a discussion of subtleties, see 
[[3~2"|]. But in experiments there are very large b —> c backgrounds and we can therefore 
only measure the portion of the rate near the lepton momentum endpoint; that is, only 
above some p min . The b — > X c £v,X u £v spectra are shown in the left half of Figure |]. The 
difficult calculation is the fraction of the leptons above a certain momentum cut. No only 
do we rely on local duality now, but terms of order 1 / (Mb — 2p e min ) can enter, spoiling 
convergence. In more physical terms, we require the detailed shape and normalization 
of the spectrum near the endpoint. 

Theory can profitably relate the endpoint b —* u£v rate to the observed b — » sy 
spectrum, since they are smeared by a common non-perturbative structure function 
[ [3~3"| , [3"4"| ], up to corrections of order Agc£>/ m fo (see p5| ] for a review). We can extract the 



structure function from b — > sy and then use this to predict the fraction of b — > u£\ rate 
above the experimental lepton momentum cut. There is still an active debate concerning 
the details of the particular methodology we employ [ ^6| ] . 

A neural net is used for continuum suppression and the signal region in lepton mo- 
mentum comprises 2.2 < pi < 2.6 GeV/c. We have lowered our cut from 2.3 GeV/c to 
increase the rate. The data are shown in the right half of Figure ^ We observe good sub- 
traction for pi > 2.6 GeV/c, and extract ( 1 874 ± 1 23 ± 326) B -> X u £v events. This yields 
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FIGURE 3. Left: typical prediction of the lepton spectra for b — > c^v and b — > m£v (note the xlO here!). 
Right: inclusive leptons in data near the endpoint. Plot a) gives the raw spectrum and shows the continuum 
(shaded) and b — » c (open histogram) contributions; plot b) shows the extracted efficiency-corrected 
b — » wfv rate. 



a partial branching ratio (before radiative corrections) of A® ;( /,(2.2 — 2.6GeV/c) = 
(2.35 ±0.15 ±0.45) x 10~ 4 . Systematics include variations of form factors and heavy 
charm states in Monte-Carlo modeling of b — > c£v backgrounds. 

To extract \V u b | , we start with an expression for the inclusive rate derived in the upsilon 
expansion [j3~7|]: 

\V ub \ = [(3.06±0.08±0.08)xl0- 3 ]x[(^/0.001)-(1.6ps/x s )] 1 /2 

The required <B u b = B(B — > X M ^V) is related to the observed rate A®^ in momentum 
window (p) by A$ u b(p) = F u (p) (B u b. The b — > 57 spectrum will provide our prediction 
for F u (p), which is simply a properly normalized integral of the observed portion of 
the solid curve on the left of Figure [|. Using b — > sy data with 1.5 < Ey < 2.8 GeV, 
we fit the shape function [0] to various parameterizations. We use this to determine 



[[38]] F u (2.2 - 2.6 GeV/c) = 0.138 ±0.034. Consistent results are obtained by the more 
model-dependent method of fitting the spectrum to parameters in the Ali-Greub specta- 
tor model [|8p and feeding this information into the ACCMM model [|33|, [TOp. 
Our preliminary result is: 

IV^I = (4.09 ±0.14 ±0.66) x 10~ 3 

This compares favorably with CLEO's exclusive (%/p/(0)£v analyses [ [14] , |4T[ ]: 

\V ub \ = (3.25 ±0.14+g;2J± 0.55) x 10" 3 



ZERO-RECOIL POINT OF B — > D*£v 



The newer, more inclusive methods above may have reduced model dependence in some 
sense, but quark-hadron duality is always involved. We now turn to a more traditional 
exclusive extraction of \V c b\ from B — > Z)*£v decays. 

This analysis measures the absolute rate as a function of q 2 . Both D* + £v and D*°£v 
modes are reconstructed, with 0.8 < p e < 2.4 GeV/c and 1.4 < < 2.4 GeV/c. One 
usually replaces q 2 with w = vg ■ vb, the product of 5 and D* meson four- velocities. 
This new variable is just a particular linear transform, w = a — bq 2 . HQET simplifies the 
analysis by relating the three form factors present to one universal Isgur-Wise function. 
In fact for B — > D*£\, the point w = 1 (corresponding to maximum q 2 ) has no O (1/M) 
corrections [fi2||. Physically, this occurs when vg = V£>, that is at zero-recoil where the 
D* is at rest relative to the B. This is the favorable place for a precise extraction of \V c b\. 

Background discrimination is accomplished by examining the angle, Bb-d*£, between 
the B and the D*£ system. This variable is similar to the familiar missing-mass. Since 
cosQ B _ D *£ is calculated from four vectors it can be unphysical and its shape can be 
helpful in disentangling the many types of background. Determining the rate for each 
bin in w requires a detailed fit like the examples on the left of Figure The right 
half of the figure shows both the fit rate and, after efficiency and kinematic factors, the 
extracted form-factor as a function of w. Inclusion of the D*° mode adds significantly 
more efficiency at zero recoil as compared to using D* + alone. 

We determine the preliminary branching ratios: 

<B(B°^D* + £v) = (5.82 ±0.19 ±0.37)% 
<B{B'^D* Q £v) = (6.21±0.20±0.40)% 

The intercept at zero-recoil and slope of the HQET form-factor (with the curvature 



related to the slope by dispersion relations [ f43j [44]]) are preliminarily determined as: 



F(l)|Vd,| = (4.22±0.13±0.18) x 10~ 2 
p 2 = 1.61±0.09±0.21 



Significant systematics include efficiency (especially for slow pions), the intermediate 
branching ratios, backgrounds, and (especially for p 2 ) the ratios of the three D* form 



factors Ri and R 2 Using F(l) = 0.913 ±0.042 we extract: 
\Vd,\ = (4.62±0.14±0.20±0.21) x 10~ 2 

where the errors are statistical, systematic, and theoretical, for a net 7% precision. 

A comparison of recent results for the slope and intercept of the B — > D*£v form factor 
is shown in Figure ^. The differing shape of the CLEO ellipse is due to an interaction of 
the lepton momentum cut with variations of the form-factor ratios R\ and R 2 within their 
errors. There is a disagreement at about the two sigma level. One difference in technique 
involves the D*X£v background; CLEO includes these in the cos0g_£>*^ fits to the data, 
while LEP analyses use a model constrained to other LEP results on B — > D*X£v. 
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FIGURE 4. Left: Signal and backgrounds for B -> D*tv in the bin 1.10 < w < 1.15 for both D*+ 
and D*° modes. Note that backgrounds can be constrained from rates in the non-physical region of the 
cos Qg^o' t- Combinatoric background refers to fake D*; (un)correlated background refers to cases where 
the real D* and I are (not) from the same B decay. Right: The top two panels show the fit yields in each w 
bin for each D* mode, while the lower panel shows the extracted form factor with the final fit overlaid. 
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FIGURE 5. A graphical compilation of the present result along with published results from LEP 
experiments p7|]. 'OPAL inc' denotes an analysis using partial reconstruction of D* + £v. Ellipses indicate 
A% 2 = 1, including systematic errors. 



THE FUTURE 



There are several related analyses in progress at CLEO. These include B — > X c £v lepton 
spectrum moments (< Ei >) which will provide another band in the A — "k\ plane. Both 
low-background tagged (di-lepton) and higher statistics untagged analyses are being 
pursued. We also have more statistics for the exclusive B — ► D*£v analysis. 

Branching ratios and form-factor investigations for the B — > (%/p/(a)£v modes used 
for \V u b\ are underway using neutrino reconstruction. We will further address inclusive 
measures of \V u b\ with inclusive leptons, but now making full use of kinematics. Instead 
of singling out the lepton momentum, one can use quantities such as q 2 and the recoil 
mass [|48|]. We hope to accept a larger portion of rate while controlling background, thus 
reducing uncertainty on fraction of the b — > u£v rate observed. 

Results for \V c b\ from D* + e + \ and from <B (B — > Xe~v) [|5(J have recently been 
presented by Belle. Future results from both Belle and BaBar will be of great interest. 



CONCLUSION 

We have measured the photon spectrum from b — ► sy decays and the hadronic mass 
moments from B — > X c £v. Using HQET and $ (B — > X c £\), we extract \V c b\ with more 
controlled theoretical systematics, but still subject to duality issues. New techniques 
relating studies of the lepton endpoint using a structure function constrained to b — > sy 
photon spectrum allow us to extract \V u b\. Lastly, we update our \V c b\ result using the 
B — > D*£v rate extrapolated to zero recoil. Ongoing analyses will extend all of this work 
further. Our emphasis is on using a variety of techniques, with differing systematics. 
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